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We study the effects of interspecific correlations in a biological coevoiution model in which or- 
ganisms are represented by genomes of bitstrings. We present preliminary results for this model, 
indicating that these correlations do not significantly affect the statistical behavior of the system. 



I. INTRODUCTION 



The dynamics of biological coevoiution poses many problems of interest to the statistical-physics and complex- 
systems communities El 13 ■ Recently, we studied a coevoiution model that is a simplified version of the one introduced 
by HaU, et al. HSHlg. Our model, in which individuals give birth, mutate and die, displays punctuated equilibria- 
like, quiet periods interrupted by bursts of mass extinctions Interactions in this model were given by a random 
interaction matrix. Here we report on a modified version of the model, in which we have added correlations to the 
interaction matrix in order to increase the biological realism. The modified model is compared with the original one 
to assess the effects of correlations. 



II. MODEL 

We used a bitstring genome of length L to model organisms in the Monte Carlo (MC) simulations In this 

sense, each genotype, which is just an L-bit number, is considered a different haploid species. Therefore, the terms 
"genotype" and "species" are used in the same sense in this paper. We denote the population of species / at a discrete 
time t as nj{t), and several of the 2^ possible species can be present at the same time in our "ecosystem." All species 
reproduce asexually (by cloning) in discrete, non-overlapping generations. In each generation t, every individual of 
species / is allowed to give birth to F offspring with a probability Pj. Whether it reproduces or not, the individual dies 
at the end of the generation, so that only offspring can survive to the next generation. The reproduction probability 
for an individual of species I is given by H, 0, IE ly| 

Pii{nj{t)}) = 1 + exp [~ E,7 Mjjnj{t)/Ntot{t) + N,ot{t)/No] ' 

The Verhulst factor A'o represents the carrying capacity of the "ecosystem" and prevents the total population TVtot (t) = 
'^jni{t) from diverging to infinity jllj. M is the interaction matrix, in which a matrix element Mij represents the 
effect of the population density of species J on species /. A positive Mjj means that species / benefits from species 
J, while a negative Mjj corresponds to a situation in which species / is harmed or inhibited by the presence of J. 
The form of the interaction matrix M is discussed in the next section. 

In each generation, all individuals undergo mutation with a probability /i. If an individual is chosen to mutate, one 
bit in its genome is picked randomly and flipped. Since we consider different genotypes as different species, mutations 
lead to speciation, i.e., creation of another species. 
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FIG. 1: Correlation functions and matrix-element distributions after averaging, (a) Theoretical (solid) and numerical (dashed) 
results for the correlation functions for 2^^ x 2^^^ interaction matrices of the single-correlated (lower) and double-correlated 
(upper) models. The theoretical and numerical results agree, (b) Distributions of matrix elements: uniform and Gaussian 
uncorrelated (solid curves with square and Gaussian distributions, respectively), single-correlated (dashed curve) and double- 
correlated (dot-dashed curve) models with L = 8. The approximately Gaussian distributions for the two correlated models and 
the distribution for the Gaussian-uncorrelated model practically overlap 



III. THE INTERACTION MATRIX 



In this study, the interaction matrix M was set up in four different ways. In the first case, all off-diagonal elements 
were randomly and uniformly distributed on [—1, 1], while all diagonal elements were set to zero. This corresponds to 
the case in Refs. jHSl- We shall call this the uniform-uncorrelated model. In the second case, we added correlations 
between interaction constants of similar species to make this model more realistic. In a real ecosystem, two different but 
closely related species X and Y interact with another species Z in a similar way. Therefore, the interaction constant 
Mxz should be positively correlated with Myz- To implement these correlations, we modified the interaction matrix 
by averaging all terms over their nearest neighbors (nn) in Hamming space. Thus, 
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where Mjj are independent variables uniformly distributed on [—1, 1], and nn{I, J) are those bitstring pairs that differ 
from the bitstring pair (/, J) by one bit (a Hamming distance of one). A square-root appears in the denominator 
because we multiply the average with the square-root of the normalization constant in order not to change the standard 
deviation of the matrix elements. 

To investigate the results of longer-range correlations we also modified the random interaction matrix by averaging 
all terms over their nearest and next-nearest neighbors (nnn) in Hamming space: 
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where nnn{I, J) is the set of bitstring pairs that differ from (/, J) by two bits. We shall call the nn-averaged and 
nnn-averaged models single-correlated and double-correlated, respectively. 

As a result of the central limit theorem, after averaging, the distributions of the matrix elements in both the single 
and double-correlated models take an approximately Gaussian form with the same standard deviation as the uniform 
distribution of the elements of the initial, random matrix (see Fig. E^)- In order to see whether a possible difference 
in the results for different models is due to the correlations or the matrix-element distributions, we also set up another 
uncorrelated interaction matrix. In this fourth model, the uncorrelated, random matrix elements are distributed with 
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Lifetime (generations) 

FIG. 2: Normalized histograms of species-lifetimes based on simulations of 2^^ generations each: uniform-uncorrelated (solid 
curve), Gaussian-uncorrelated (short-dashed curve), single-correlated (dot-dashed curve) and double correlated (long-dashed 
curve). The distributions for the models with a Gaussian matrix-element distribution overlap (within the margin of error), and 
they differ significantly from the lifetime distribution for the uniform-uncorrelated model. The histograms exhibit a power-law 
like decay with an exponent near —2. Results are averaged over eight runs each 

a Gaussian distribution with the same standard deviation, -^/l/S, as in the correlated models. We shall call this the 
Gaussian-uncorrelated model. 

The correlation functions and the distributions of interaction-matrix elements for all the models are shown in 
Fig. n] The steps in the correlation functions are a result of the chosen metric. We use a city-block metric in which 
the Hamming distance between two matrix elements M/j and Mkl is given by H{I, K) -\- H{J, L) where H{I, K) is 
the Hamming distance between two i-bit bitstrings, / and K. 

IV. SIMULATION RESULTS 

The interaction matrix is set up at the beginning and is not modified during the course of the simulation. For all 
models, we performed eight sets of simulations for 2^^=33 554432 generations with the same parameters as in Refs. 
1^0: genome length L = 13, mutation rate ^ = 10^"^ per individual per generation, carrying capacity A^o — 2000, 
and fecundity F = 4. We began the simulations with 200 individuals of genotype 0. In order to compare the models, 
we constructed histograms corresponding to the species- lifetime distributions. The lifetime of a species is defined as 
the number of generations between its creation and extinction. As seen in Fig|21 the lifetime distributions of the 
correlated and Gaussian-uncorrelated models overlap within the margin of error, and they differ significantly from the 
lifetime distribution of the uniform-uncorrelated model. They all exhibit a power-law like decay with an exponent near 
—2. The correlations between the matrix elements do not seem to affect the behavior of the lifetime distribution to a 
statistically significant degree, at least not for the relatively weak correlations that were introduced here. On the other 
hand, changes in the marginal probability density of the individual matrix elements do have statistically significant 
effects, even though gross features, such as the approximate 1/x^ behavior of the species-lifetime distribution are not 
changed. Similar conclusions are reached also for other quantities that we studied. In particular, the power-spectral 
density of the Shannon- Wiener species diversity index shows 1// noise [^|^ with an overall intensity that depends 
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more on the marginal matrix- element distribution than on correlations in M. 

Although we have tested only weak correlations, and so our conclusion is only preliminary, there appears to be no 
disadvantage in using a random interaction matrix to model such an "ecosystem." This has obvious computational 
advantages as it makes it possible to simulate systems with larger numbers of completely different species. 
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